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ABSTRACT 

We analyze the reliability of oxygen abundances and ionization parameters obtained from 
different diagnostic diagrams. For this, we compiled from the literature observational emis- 
sion line intensities and oxygen abundance of 446 star-forming regions whose O/H abundance 
was determined by direct estimation of electron temperature. The abundances compiled were 
compared with the values calculated in this work using different diagnostic diagrams in com- 
bination with results from a grid of photoionization models. We found that the [Om]/[On] 
vs. [Nn]/[On], [Oml/HjS vs. [Nii]/[Oii], and ([Oin]/H/3)/([Nn]/Ha) vs. [S n]/[S in] diagnos- 
tic diagrams give O/H values close to the r e -method, with differences of about 0.04 dex and 
dispersion of about 0.3 dex. Similar results were obtained by detailed models but with a dis- 
persion of 0.08 dex. The origin of the dispersion found in the use of diagnostic diagrams is 
probably due to differences between the real N/O-O/H relation of the sample and the one as- 
sumed in the models. This is confirmed by the use of detailed models that do not have a fixed 
N/O-O/H relation. We found no correlation between ionization parameter and the metallicity 
for the objects of our sample. We conclude that the combination of two line ratio predicted by 
photoionization models, one sensitive to the metallicity and another sensitive to the ionization 
parameter, which takes into account the physical conditions of star-forming regions, gives 
O/H estimates close to the values derived using direct detections of electron temperatures. 
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1 INTRODUCTION 

Oxygen abundance estimates in star-forming regions play a cru- 
cial role in the understanding of galaxy evolution. For example, 
oxygen radial gradie nts in spiral galaxies obtained by Hn re 



gion observations (e. g.[Stanghellini et al. l20ld : Kewl ev et al 
Bresolin et ail 120091: iKrabbe et alj 120081 : iDors & Copettil 



i Bres 
kem 



2010; 



2005 



Kennicutt et alj|2003l) are ess ential to test chemical evolution mod- 
els (see lMolla & Dfaal2005l) a nd to investigate the effect of envi- 
ronment on galaxy interacti ons fellison et afbOlollDors & Copettil 
2006; Slril lmanetal.il 19961) as well as the mass-metallicity relation 
of ga laxies (e.g. Pilvugin et alj |2004| ; IPerez-Montero & Continil 
2009). Likewise, oxygen abundance estimates in metal-poor galax- 
ies are also important to test theories of chemical evolution of 
;alaxies because these are the least chemically evolved objects 
Kunth & Sargentl 19831) . 

Unfortunately, for the most of star-forming regions, only 
collisionally excited emission-lines (CELs) in the optical are 
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bright enough to be used for the derivation of elemental abun- 
dance. CELs are temperature sensitive, thus, only an accurate 
determination of the metallicity can be achieved from the previous 
estimation of the electron temperature (this method will be 
called r e -method) using, for instance, the ratio of different CELs 
[Om](/l4959 + ,15007)/,i4363 which are weak or unobservable 
in st ar-forming regions with high meta llicity and/or low excita- 
tion toors et al.1 120081 iDfaz et alj|2007l) . In these cases, oxygen 
abundances can be obtained by empirical (i.e. using oxygen deter- 
minations via r c -method) or theoretical (i.e. using photoionization 
models) calibrations between oxygen abundances and more easily 
measured line ratios (hereafter strong-line methods). The oxygen 
abundance i ndicator flp=([On ]/i3727+[0 m]A4959,A5001)/Uj3 
proposed by IPagel et all d 19791) has found large acceptance 
in this context and several authors have calibrated th is line 
ratio with O/H abundance (e ,g. lEdmunds & Pagel 



iDopita & Evansl 1 19861 : IPilvuginl l200ll : IDors &C opetti 



1984 



2005 



among others). Additional O/H indicators based on other emission 
lines such as A' 2 = [Nn]^6584/Ha 
1 19941) , [Nn]^6584/[Om]/l5007 



(Storchi-Ber 


gmann ct al. 


jAlloin et al. 


19791). 



Dors et al. 



5 23 =([S n]^6716,6 7 31+[S n],U9069,9532)/Hff 
dVflchez & Estebanl 1 19961 : iDfaz & Perez-Monterol |200Ch have 
also been suggested (see also lKewlev & Dopitdl2002h . However, 
distinct methods or distinct calibrations of a same oxygen indicator 
provide different oxygen values with discrepancies up to 1.0 dex 
dKewlev & Ellisordl2008l : iRupke et all 20081 : iDors & Copettl 20051 : 
Kenni cutt et al 20031 ). Currently, the large number of direct oxy- 
gen estimates available in the lit erature have help ed to investigate 
this discrepancy. For example, I Yin et all {2007) determined the 
gas-phase oxygen abundance for a sample of 695 galaxies and H n 
regions using the r c -method and compared these determinations 
with the ones via R 23 , N 2 , ([Nn]/l6584/Ha)/([Om]15007/HjS), 
and [S ii](/16717+/16731)/Hq'. They found that among the indices 
above, the N 2 provides more consistent O/H abundances when 
compare d with the ones via 7" e -method . Similar analysis wa s also 
done by IPerez-Montero & Dfazl d2005l) . Iliang et all J2006h . and 
iNagao et alJfcOOeh . 

The studies above analyzed strong-line methods based mainly 
on one line ratio, such as the R 2 i, N 2 , among ot hers. In princi- 
ple, th e use of diagnostic diagrams, suggested by iBaldwin et al. I 
(1981) to separate objects according to their primary excitation 
mechanisms, containing line ratios strongly dependent on the de- 
gree of ionization and on the metallicity of star forming regions, 
can improve the accuracy of strong-line methods. Although a large 
number of these diagrams have been applied to estimate oxy- 
gen abundances and ioniz ation parameters of star forming regions 



[OIII]/Hp 



(e.g. Levesque et al. 2010, Viironen et al. 2007, Kewley et al. 2001, 
iMcGaugh et. al II 199lL IDors et alj|2008l |Popita & Evansl 1 198rJT a 
comparison of oxygen estimates obtained from these diagrams and 
r c -method is unavailable in the literature. 

Another important issue related with the determination of 
metallicity using strong-line methods is the relation between the 
ionization parameter and the meta l licity, which is still contro- 
versial. For example, Gar nett et al.l d 1997b . in a study of the in- 
terstellar abundance gradient in NGC 2403 found that any cor- 
relation between ionization parameter and abundance must be 
a weak one. This result is in agreement with that obtained by 
iKennicutt & Garnettl d 19961) from measurements of the [S II]/[S 
III] ratio in 41 H II regions in M101 and with the on e foun d 
for three barred g alaxies established by IDors & Copetti I d2005h . 
On the other hand, iBresolin et al.l d 19991) found that in metal poor 
disk Hn regions the ionization parameter is about 4 times larger 
than in H 11 regions w i th sol ar m etallicity. Thi s relati on was also 
found by iMaier et al] d2006t) and iNagao et all Ij2006) for a sam- 
ple of galaxies. This subject is important to study the association 
of the mass-metallicity relation with the mass-age relation in local 
gal axies or the relation between gas metallicity and stellar metallic- 
ity dNagao et al.ll2006h . Additional analysis using several methods 
would help to solve this disagreement. 

In this paper, we employed a grid of photoionization models 
and data compiled from the literature in order to estimate the oxy- 
gen abundance using diagnostic diagrams to compare with those 
via re-method, as well as to investigate the relation of the ioniza- 
tion parameter with the metallicity. Detailed photoionization mod- 
els are also built to produce more precise determinations of these 
parameters. In Section [2] we describe the observational data used 
in the analysis. The modeling procedures are presented in Sect.[3] A 
description of the diagnostic diagrams employed is given in Sect.|4] 
The results and discussion are presented in Sects. |5]and|6] respec- 
tively. A conclusion of the outcome is given in Sect. [7] 
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Figure 1. Comparison of some predicted line ratio intensities from pho- 
toionization models with nebular and stellar metallicity matched (red) and 
no matched (black). 



Table 1. Emission line ratios considered 



Symbol 


Definition 


R23 


([O n]i3727+[0 mM4959,i5007)/H/3 


[O mynp 


[O HI]A5007/He 


[Oiii]/[Oii] 


[O my5007/[On]/l3727 


[NIIl/Ho- 


[NII]i6584/Hff 


[NII]/[OII] 


[NII]i6584/[OIIW3727 


[S II]Ha 


([S II]i6716+i6731)/Ha 


[SUMS III] 


([S II]i6716+i6731)/([S III] ^9069+^9532) 



2 OBSERVATIONAL DATA 

Observational emission line intensities of a sample of H n galaxies 
and H 11 regions and the oxygen abundance computed using the T e - 
method were compiled from literature. 

The emission lines considered in our analysis are listed 
in Table [T] Thi s compilation inclu d es data of Hn re - 
gions obtaine d by Bresolin et al. (2009), Bresolin et al. (2005 
IBresolin et ail d2004l). IBresolin! OOOTTTlLee & Skillmanl i200- 
and IKennicutt et alj 120031). T he data on Hn galaxies were ob- 
tained by llzotov et alj d2006l), iLee et alj (Eoolf). Ilzotov & Thuanl 



JL 

d2004l), IVflchez & Iglesias-Paramol d2003l) , lHagele et alj d2008E 
and iGuseva et alj d2000l) ~ The sample consists of 446 objects (86 
H 11 regions and 360 H 11 galaxies) whose O/H abundance are in the 
range 7.0 < 12 + log(0/H)r c < 9.0 and represents practically the 
entire metallicity range of star-forming regions (see IPilvugin et al.l 
l2004al) . The objects of the sample have z < 0.07 and measurements 
corrected by dust extinction and no AGN and gas shock contribu- 
tions are present in their ionization. 



3 PHOTOIONIZATION MODELS 
3.1 Model Grid 

To enable the estimation of the oxygen abundance and the ioniza- 
tion parameter using diagnostic diagrams, we employed the pho- 
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Table 2. Final parameters of the detailed photoionization models for H n regions observed bv lKennicutt et aljfe003l) 



H ii region 12+log(0/H) log(N/0) log(S/0) log U (cm" 3 ) Age(Myr) 



H1013 


8.48 


-0.72 


-1.34 


-1.94 


47 


2.5 


HI 105 


8.71 


-1.06 


-1.64 


-2.35 


248 


1.0 


H1159 


8.80 


-1.28 


-1.94 


-2.55 


5 


1.0 


HI 170 


8.07 


-0.83 


-1.31 


-2.60 


9 


2.5 


HI 176 


8.16 


-0.76 


-1.38 


-2.05 


33 


2.5 


H1216 


8.00 


-1.28 


-1.64 


-2.52 


33 


1.0 


H336 


8.75 


-0.80 


-1.58 


-2.68 


15 


2.5 


H409 


8.53 


-1.15 


-1.58 


-2.37 


213 


1.0 


H67 


8.00 


-1.15 


-1.64 


-2.82 


10 


2.5 


N5471-D 


8.10 


-1.15 


-1.64 


-2.4 


110 


2.5 


SDH323 


7.76 


-1.45 


-1.61 


-2.91 


61 


2.5 



Table 3. Observed intensity lines and those predicted by our detailed models 



H ii region 


[OII]<13727 


[OIII]/14363 


[OIII]i5007 


[NII]i6584 


[S II]/16720 


[S III]i9069 + ,19532 




Obs. Mod. 


Obs. Mod. 


Obs. Mod. 


Obs. Mod. 


Obs. Mod. 


Obs. Mod. 



H 1013 


188±10 


185 




0.20 


103±5 


105 


64.6±3.4 


67 


28.8+1.1 


26 


131.6±7.0 


171 


H1105 


185±10 


188 


1.4+0.1 


0.74 


316+17 


317 


33.4+1.8 


34 


24.1 ±0.9 


27 


126.6±11 


139 


H1159 


198±10 


194 


1.9±0.4 


0.63 


317+17 


316 


23.6+1.3 


25 


29.8±1.1 


27 


97.0 ±5.5 


66 


H1170 


308±16 


295 


1.6±0.2 


1.67 


201+11 


194 


44.0±2.3 


41 


56.7±2.1 


51 


170.0±9.3 


136 


H1176 


160±8 


153 


2.4±0.3 


2.74 


369±20 


357 


21.2+1.1 


23 


23.0±0.8 


19 


1 13.5+6.1 


137 


H1216 


151±8 


149 


4.7±0.3 


6.67 


473±25 


478 


11.0±0.6 


9 


18.9±0.7 


18 


83.0±4.6 


63 


H336 


178±9 


183 




0.02 


23+1 


26 


95.9+5.1 


99 


56.8±2.1 


58 


107.0±5.7 


102 


H409 


218±12 


212 


2.3±0.2 


1.39 


370±20 


359 


27.3±1.4 


24 


31.2+1.1 


28 


90.1 ±4.9 


141 


H67 


244±13 


248 


3.5±0.5 


4.67 


342±18 


356 


16.3±0.9 


19 


26.3±1.3 


29 


92.1±5.7 


62 


N5471-D 


137±7 


140 


8.0±0.4 


7.10 


578+31 


574 


8.5±0.5 


11 


20.6±0.8 


18 


75.7±4.3 


75 


SDH323 


194±10 


198 


5.5±0.9 


3.7 


227±12 


234 


7.9±0.7 


7 


20.8±1.4 


23 




42 



toionization code Cloudy 8.00 dFerlandl2002l) to build a grid of 
models covering a large space of nebular parameters. In these mod- 
els, a stellar cluster was assumed to be responsible for the ionization 
of the nebulae, with a spectral energy d istribution (SED) obtained 
using S tarburst99 Jlxitherer et al.|[l999h . We built models with 
stellar cluster formed by instantaneous burst with Salpeter initial 
mass function (a = -2.35), lower and upper stellar mass limits of 
0.1 M G and 100 M Q , respectively, and age of 2.5 Myr. Other papers 
that have considered stellar clusters as ionizing s ources in order 
to reproduce strong forbidden lines of H n reg ion ( Dors & Copetti 
al.lbOCKtr 



2005;lDopita et 



gion ( L) 

IStasinska & Izoto v 2003; Bresolin et al 



1999; Copetti et alj 119851) have derived about the same age for 



star-forming regions (i.e. 1-3 Myr). Similar ages have also been 
found from optical ph otometric data of giant Hn regions (e.g. 
iMawa & Prabhu 1 1 1996b . Selection effects may explain this lim- 
ited range of ages. Hn regions younger than about 1 Myr are 
difficult to be detected in the optical, because they are gener- 
ally embedded in dusty molecular clouds which cause consider- 
able optical extinction. Nebulae older than about 5 Myr are also 
difficult to be obse rved because their original massive stars have 
cooled or are dead jDopita et al. Il200oi iGarcfa- Vargas et ail 19961 ; 
ICopettietalJll985h . We used the stellar evolution models from the 
Genev a group with high m ass-loss rates and without stellar ro- 
tation dMevnetetal.l[l 994). The non-LTE atmosphere model of 
IPauldrach et alj 1 2001 ) was assumed in the models. If LTE atmo- 
sphere model is assumed instead of non-LTE model a lower ioniza- 



tion degree is produced in the hypo thetical nebulae jDors & Copettil 
120031 ; IStasinska & Schaererlll997h . This would affect mainly ion- 
ization parameter rather than metallicity determinations via strong- 
line methods. The models were built having ionization parame- 
ter ranging from log U = -1.5 to -3.5 (with a bin size of 0.5 
dex), metallicities (traced by the oxygen abundance) Z= 0.04, 0.02, 
0.008, 0.004, 0.001, plane-parallel geometry, and electron density 
of N e = 200 cm" 3 . This electron den sity value is typical of not 
evolved H 11 regions dCopetti et al. 2000). 

The abundances of heavy metals in the nebula is scaled lin- 
early to the solar metal composition through the comparison of 
the oxygen abundances, with the exception of the N abundance, 
w hich was taken from the r e lation log(N/O)=log(0.034+ 1 20 O/H) 
of IVila-Costas & Edmunds! jl993h. T he so lar composition (Z= 
0.02) refers to lAllende Prieto et all d200ll) and correspond to 
12+log(0/H)= 8.69. The presence of internal d ust was considered 
and the grain abundances dvan Ho of et al. 200]]) were also linearly 
scaled with the oxygen abundance. To take depletion of refractary 
elements onto dust grains into account the abundances of the ele- 
ments Mg, Al, Ca, Fe, Ni, and Na were red uced by a factor of 10, 
and Si by a factor of 2 dGarnett et al .11 1995b relative to the adopted 
abundances in each model. 

The solar metallicity for the stars from the Geneva evolution- 
ary tracks, which correspon ds to the old solar oxygen abundance 
value [12+log(0/H)= 8.87 dGrevesse & Sauval II 19981) 1 . is higher 
than the value adopted for the nebular component. This produces 
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-1 -1.5 -1 -0.5 



log([NII]/[OII]) log([SII]/[SIII]) 

Figure 2. Diagnostic diagrams containing observational data taken from literature (see Sect. [2} and results of the grid of photoionization models (see Sect[5}- 
The solid lines connect curves of iso-Z, while the dotted lines connect curves of iso-f/. The values of log U and Z are indicated. Squares are the observational 
data. The line ratio used in each plot is defined in Table[TJ The typical error bar (not shown) of the emission line ratio is about 10 %. 



an imperfect match between gas and st ar metallicity i n our pho- 
toionization models. As pointed out by iDopita e t al. ( 200^), the 
main effect of this is that the computed stellar UV photon field 
result slightly softer. To investigate how this discrepancy affects 
our results, we built a model with a perfect match between neb- 
ular and stellar metallicities and compared the result with another 
model whose metallicities are in disagreement. This latter was built 
using a SED with 12+log(0/H)=8.69 obtained by linear interpola- 
tion of the spectra with Z= 0.02 and Z= 0.008. In Figure [T] we 
show the histogram with the comparison of some emission line 



intensities ratio predicted by these models. We can see that, with 
exception of [OIII]/[OII] and [OIII]/H/3, the intensities of the 
majority of the line ratios show little variation when the stellar 
metallicity atmosphere change. Similarly, the predicted value of 
the R23 index (not shown in Fig. [TJ shows a small deviation, i.e. 
about 12%, which corresponds to variations in the oxygen abun- 
dance from calibrations using this line ratio by only 0.02 dex. 
Thus, the disagrement between the nebular and stellar metallicity 
have little influence on the Z determinations from strong-line meth- 
ods and can only affect U estimates in diagnostic diagrams which 



Analyzing derived metallicities and ionization parameters from model-based determinations 5 



use th e [OIII]/[OII] and [OIII]/Hj8. Interestingly. iDopita et all 
d2006l) found that critical line ratios changed by 0.1 dex or less, 
except for the [O i]/l6300/Ha ratio, when a test model with a 0.4 
Zq spectral synthesis cluster model from Starburst99 embedded 
in a nebula with 1.0 Zq (12+log(0/H)= 8.66) is run. Along the pa- 
per the solar abundance ado pted refers to 12+log(0/H)= 8.69 from 

lAllende Prieto et al.1 J200lh . 

These models are similar to the ones of Dors & Copetti d2006l) 
and have b een successful in describing observational data of H n re - 
gions (see lDors et al .l2008l ; lKrabbeet alj200 8. Krabbe et alj2007h . 



3.2 Detailed models 

In general, grids of photoionization models are built assuming 
a fixed N/O-O/H relation. However, this constancy can yield 
large uncertainties in O/H est imates via strong-line methods 
( Perez -Montero & Continil l2009h . This problem can be circum- 
vented by the use of detailed photoionization models. To analyze 
the source of these uncertainties, we built detailed models in order 
to reproduce the observational emission line intensities of 1 1 H n 
regions (see Table[2} located along the disk of the galaxy M 101, 
observed by iKennicutt et alj ( 120031 ), and we compared our esti- 
mates with O/H and U values from other methods. These objects 
were selected because they cover the wide range in metallicity and 
ionization parameter considered in this paper. We computed indi- 
vidual models for each object adopting the following methodology. 
Firstly, a model for each region was built by initially guessing the Z 
and U values derived from a comparison between the grid of pho- 
toionization models shown in the diagnostic diagram [Om]/[On] 
vs. [Nii]/[Oii] (see Figure|2j and the observational data. The elec- 
tron density of each model was considered to be that computed uti- 
lizing the task temden of the package IRAF, where we consider the 
sulfur ratio [S n]/16716 /[S n]/l6731 and elec tron temperature for the 
+ ion measured bv lKennicutt et alj (2003). The stellar cluster was 
assumed to have an age of 2.5 Myr, M up = 100 M Q and metallic- 
ity was matched with the closest one nebular assumed in the mod- 
els. Then, we ran new models ranging the O/H and U values by 
0.3 and 0.5 dex, respectively, with a step of 0.1 dex. From this se- 
ries of models we selected a model which produced the smallest 

TiXf = xfoil]/H/3 + ^foiII]/H£' w here Xi = U' bs. ~ ^pred.)~/'obs.' -'obs. 

and /p red are the observational and predicted intensities of the line 
ratios, respectively. Another series of models was computed con- 
sidering the O/H and U values found by the criterion above but 
ranging the N/H and S/H abundances by 0.3 dex in order to re- 
produce the intensities of the [Nn]/i6584 and [Sn]/16720 emission 
lines. The satisfactory solution is found when / prc d. reproduces / D b s . 
within the observational uncertainties and the model has the small- 
est TjX] ~ x\ouyup + xlomyap + Xymiyup + ^fsii]/H/r ^ n some cases 
no satisfactory solution was reached considering the age of the ion- 
izing cluster of 2.5 Myr. For these, it was necessary to assume an 
age of 1 Myr because their observed emission lines could only be 
reproduced by means of a harder spectral energy distribution. In 
Table|2] we present the final parameter obtained for the models. 



4 DIAGNOSTIC DIAGRAMS 

We employ six diagnostic diagrams containing predicted and ob- 
served emission line ratios sensitive to Z and U. The diagrams con- 
sidered are described below. 

• [Om]/[On] vs. [Nii]/[Oii] — Diagnostic diagram suggested 




6000 7000 

Wavelength (A) 



Figure 3. Ratio between observed and predicted emission line intensities of 
some Hn regions located in M 101 observed bv lKennicutt et al . 12003). 



by IDopita et al. I d2000l) . where the [O ra]/[0 n] has a strong de- 
pend ence on U, on the effective temperature of the ioniz ing stars 
(e.g. iDors & Copettil 120031: iPerez-Montero & Diazll2005l) a nd on 
the metallicity dKewlev & Dopitall2002l ; IDopita et al. Il200d) . The 
[Nn]/[On] correlates st rongly with Z above to Z > 0.04 Z 
dKewlev &popitall2002l). it is also de pendent on the N/O abun- 
dance ratio dPerez-Montero & Contim1l2009l) and it is almost inde- 
pendent of U dKewlev & Dopitall 20021) . 

« [Nn]/Hq vs. [S ii]/Hq' — This diagram was proposed by 
IViironen et al.l d2007l) to estimate t he metallicity, where both line 
ratios are dependent on U and Z dStorchi-Bergmann et al.|[l994t 
iKewlev & Dopit al [20021) . iMazzuca et al. 1 2006) pointed out that 
the use of diagnostic diagrams using [N ii]/Hq' can yield degener- 
ate values for Z, as star-forming regions with low Z and U have 
[Niij/HQ- values similar to regions with high Z and (J. In addition, 
no consistent values were found for oversolar abundances because 
the [Nii]/Hq' parameter saturates in this high-metallicity regime. 
Therefore, we did not consider models with metallicity over-solar. 
For the [S ii]/Hq', this line ratio is strongly depen dent on U and 
increa ses with the abundance for low metallicities dLevesque et al.1 

l2oToh . 

[O m]/ H 8 vs. [Nii]/[Oii] — Diagram proposed by 
d2000l). The [OmJ/HjS was suggested by 
1984) as an O/H indicator. However, due 



Dopita et al. 
Edmunds & Pagell 



to its dependence on U dDopita & Evansl [l986; McGaug h et. al I 
1 19911) a combination with another line ratio is preferable, other- 
wise, crud e O/H estimates with u ncertainties of about 0.5 dex are 
produced (Kob ulnickv et alj|l999h . 

• [N ii]/Hq; vs. [ S n]/[S m] — The [S n]/[S m] was proposed to 
be a U indicator bv lDfaz et al. Idl99ll) f or mo derate to high metal- 
licity regime (see also IDopita & Evanslll986l) . and it has little de- 
pendence on Z. The problem in using thi s line ratio is t hat it is un- 
derestimated by phot oionization models (|Garnettll989l) . especially 
for high metallicity dDors & Copettil2005l) doing any estimation in 
this regime is somewhat uncertain. 

• ( [Om]/H/3)/([Nn]/Ha) vs. [Sn]/[Sm] — IPettini & Pagell 
d2004l) showed that ([Om]/H/?)/([Nn]/Ha) is dependent on Z. Be- 
cause this line ratio is also dependent on U, we combined it with 
the [S n]/[S m] in order to minimize the uncertainties in Z determi- 
nations. 
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Figure 4. Comparison between the oxygen abundances derived using the 
re-method with those computed using the diagnostic diagrams. In each plot 
the diagnostic diagrams used to compute the oxygen abundance are indi- 
cated. The top panel of each diagram shows the difference between oxygen 
abundance via T c -method and via diagnostic diagrams. The average value 
of this difference and the dispersion are shown in each plot. Solid lines rep- 
resent the equality of the two estimates. The results for H n galaxies and H n 
regions are marked by red and black points, respectively 



• [Nii]/Hq' vs. [Oiii]/[Oii] — We investigate the combination 
between these line ratios to eliminate the problem existing with the 
use of [Sn]/[Sm]. 



5 RESULTS 

In Figure[2] the diagnostic diagrams described above containing the 
results of our grid of photoionization models and the data sample 
is shown. The majority of the observational data fall within the re- 
gions occupied by the models. However, in the diagram [Om]/H/? 
vs. [N n]/[0 n] the models predict [O m]/H/3 values lower than the 
observed ones for the low metallicity regime and high U values, 



result also found for the ([Om]/H/?)/([Nii]/Hci') line ratio. Similar 
difficulty in modeling metal poor star- forming regio n s have been 
found by other authors. For example, iDopita et al. I {2006), who 
used the same SEDs used in this paper, found that their models 
do not reproduce the o bserved emission line dia gnostic ratios of 
objects with Z < 0.4Zn. lMartm-Mani6n et all d2008h . using a com- 
bination of photoioni zation models and sets of stellar yields from 
Gavil an et al. also found that their models do not repro- 

duce iobseryational^ata_of most metal deficient H n galaxies (see 
also lFernandes et al. Il2003h . iKewlev et alj d200ll) pointed out that 
this disagreement is due to that stellar ionizing spectra are not 
hard enough in the far ultraviolet region, and inclusion of the ef- 
fects of continuum metal opacities in stellar atmospheres should 
be a way of improving the models accurac y. However, in our wor k 
we adopted stellar atmosphere models of IPauldrach et al] (1200 lh . 
which include treatments of continuum metal opacities, and still 
this disagreement can be noted. A solution for this problem seems 
to include the effects of rotation in stellar models (see discussion 
above). The number of points ranges in the diagrams because for 
some data s et no all emi s sion l ines considered were observed. For 
example, in llzotov et alj J2006l) the [O ii]/13727 is not observed in 
30% of the objects. 

For the detailed models, a comparison of the predicted and ob- 
served emission line intensities is listed in Table[3] and Fig.[3]shows 
the ratio between these. We can see that the models reproduce very 
well (with differences lower than ~ 15 %) all the observed inten- 
sities within of the observational uncertainties, with exception of 
the [S m]/19069 + /!9532 and [Om\M4363 emission lines, which are 
reproduced only for H336 and NGC5471-D; and for NGC 1170 
and NGC 1 176, respectively. Other works have also found that pho- 
toionization models are inc apable of reproduce emission line in- 
tensit ies sensitive to T e (e.g. lStasi nska & Schaerer 1997l; IOev et alj 
2000). This problem has been atributted to temperature gradients 
and/or temperature inhomogeneities in nebulae wh ich are not taken 
into account in simple photoionization models (see Stasinska 2002) 
such as the ones used in this paper. Because the [Om]/14363 emis- 
sion line has an exponential dependence on the electron tempera- 
ture, any offset between the electron temperatures in photoioniza- 
tion models and observed forbidden-line temperatures will have the 
strong effect on the reproducion of this emission line. 



5.1 Abundance determination comparison 

To check the reliability of diagnostic diagrams, in Figs. |4] and [5] 
we present a comparison of O/H obtained from these methods 
with O/H abundances via r e -method and the difference among 
these estimations. We include in Fig. [5] the results of the de- 
ta iled models for the H n regi ons in M 101 as well as the results 
of IPerez-Montero et al. who built detailed photoionization 

models in order to reproduce emission line intensities of 10 Hn 
galaxies. In each plot, the average value (D) and the dispersion (tr) 
of this difference are also presented. The O/H and U values from 
the diagnostic diagrams were obtained by linear interpolation from 
the model grid shown in Fig. [2] In a few cases double values of Z 
and/or U for a same objects are found because the models overlap 
for a given combination of these parameter. For these, the estima- 
tions were not considered in our analysis. This happened mainly 
for the [N n]/Ho' vs. [S ii]/Hq' diagnostic diagram (for about 5% of 
the points). 

The diagnostic diagrams that provide the best results are 
the [Oiii]/[Oii] vs. [Nii]/[Oii], [Om]/H/3 vs. [Nii]/[Oii], and 
([Om]/Hj3)/([Nn]/Ha) vs. [Sii]/[Sm], which give O/H estimates 
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close to the r e -method with an absolute difference of about 0.04 
dex. The lowest dispersion is found with the use of the [O iii]/H/? 
vs. [N n]/[0 n] diagram. For the majority of the diagrams, the differ- 
ence and the dispersion are larger in the regime of low metallicity 
(12+log(0/H)< 8.0). For the other diagrams this difference is about 
0.25 dex. The O/H abundances via detailed models are in conso- 
nance with the ones via T e -method for the objects analyzed and the 
dispersion derived is lower than the one obtained using diagnostic 
diagrams. 



5.2 Ionization parameter determination 



For the ionization parameter, in Fig. [6] we plotted U against the 
oxygen abundance obtained from diagnostic diagrams presented 
in Section [4] as well as those obtained from detailed models. The 
results for Hn galaxies and Hn regions are indicated by differ- 
ent symbols (red and black dots, respectively). The [O m]/H/3 vs. 
[N n]/[0 n] and [O m]/[0 n] vs. [N n]/[0 n] estimates larger U val- 
ues than the ones via other methods. There is not a clear trend of 
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U with O/H and this result is also confirmed by the detailed model 
estimation. 



6 DISCUSSION 

Comparing oxygen abundance determinations via T e -method for 
the sample of Hn regions and Hn galaxies with those based on 
strong emission lines, we found that the [O m]/[0 n] vs. [N n]/[0 n], 
[O m]/H/3 vs. [N n]/[0 ii], and ([O m]/H0)/([N n]/Hor) vs. [S n]/[S in] 



diagnostic diagrams give O/H values nearest to the r c -method, 
with differences of about 0.04 dex and a dispersion of about 0.30 
dex. This difference is about the same as the one between oxygen 
estimates via the f-method I Pilvugin 1120011) and via T e -method 
found by Lopez-Sanchez & Estebanl ( l2010h considering a sample 
of Wolf-Rayet galaxies. It is lower by about 0.15 dex than the one 
found by these authors when using only one emission line ratio 
sensitive to metallicity. 

As seen in Figs. [4] and [5] large difference are found for 
12 + log(0/H) < 8.0 (Z < 0.2 Z G ). Similar results were also found 
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by I Yin et alj d2007t) . who compared oxygen estimates via N2 and 
([O m]/H/?)/([N nJ/Ha) with those via 7Vmethod for a sample of 
695 galaxies and Hn regions. This occurs because in this regime 
of metallicity the nitrogen and oxygen have both mainly a primary 
nucleosynthesis origin, doing nitrogen emission lines to be rela- 
tively independent on oxygen abundance and consequently the use 



of metalli city indicators based on these emission lin es are not reli- 
able. (e.g. lLevesque et al.ll2010l : lDopita et al. 1120001) . 

The origin of the dispersion found by us is probable due to 
the difference between the real N/O-O/H abundance relation of 
the object sample and the one assumed in our models. In fact, 
Perez-Mont ero & Continil {2009) analyzed the dependence of N/O 
with O/H estimation obtained via the metallicity indicators using 
nitrogen line ratios and compared these estimations with the ones 
obtained via T e -method. They found approximately the same dis- 
persion as the one derived by us, and also showed that if the N/O 
ratio is taken into account in strong-line methods, the dispersion 
can be reduced by about 0.1 dex. Moreover, the scattering of N/O 
for a fixed O/H value i s larger for the low metallicity regime (see 
e.g. IPilvugin et ai]|2003l) . which introduces a larger dispersion for 
oxygen estimations in this regime, such as the one observed in our 
results. This is confirmed by the use of detailed models, for which 
the N/O-O/H relation is a free parameter, yielding a lower disper- 
sion (0.08 dex) than the ones obtained from diagnostic diagrams. 
lYin et alj {200% also obtained similar results comparing oxygen 
abundances deriv ed from r e -method and th ose via the photoioniza- 
tion models o f lCharlot & Lon ghetti 1 d200lh . 

Another important test is to verify if abundance gradi- 
ents estimates by using diagnostic diagrams agree with those 
via T c -method. For that, in Fig. [7] we show a comparison 
of oxygen gradient slope computed using the [Om]/[Oii] vs. 
[N n]/[0 11] diagram presented in Fig. [2] and those via r c -method 
for spiral galaxies M 101 M51, M33, and NGC2403 obtained 
by iKennicutt et alj f200l) iBresolin et alj J2004) . iMagrini et all 
d2007l) , and iGarnett et alTdl997t) , respectively. We can see that, 
within the uncertainties given by the linear fitting, the diagnos- 
tic diagram above yields abundance gradient consistent with the 
ones via r c -method. Again, the difference between the gradient es- 
timates is probable due to the N/O-O/H relation assumed in our 
models and the one of the galaxies. This is supported by the de- 
tailed model results, since a linear fitting on oxygen abundance 
from these, presented in Table [2] yields a gradient for M 101 of 
12+log(Q/H) = 0.90(±0.26) RJR25 + 8.77(±0.15), the same gradi- 
ents found bv lKennicutt et alj J2003t) using the Te-method. 

In general, oxygen determination obtained from strong-line 
methods, which use emission line intensities predicted by photoion- 
ization models, are overestimated up to 0.5 dex when compared 
with those obtained from T e -method jKewlev & Ellison! [2008 ; 
Dors & Copettil 120051: IKennicutt etafl |2003| ; IGarnett et alj |2004 



Stasihskal2002h . This discrepancy is attributed to the fact that pho- 
toionization codes are not realistic enough, do not treat all the rel- 
evant physical process es correctly, use inaccurate atomic data, etc 
(Kenn icuttetalj2 003). However, as seen previously, using the state 
of art of photoionization models and the combination of two line 
ratio, one sensitive to the metallicity and another sensitive to the 
ionization parameter, which does taken into account the physical 
conditions (hardness of the i onizing radiatio n and geometrical fac- 
tor) of star-forming regions dPilvugin l200ll) , minimizes the effects 
mentioned above and gives O/H estimates close to the r c -mefhod. 
As explained in Sect. [3] the match between solar abundances for 
the gas and star has little influence on metallicity indicators (i.e. 
[N n]/[0 11]) showing that the metallicity estimates by our models 
are indepent of this fact. 

The ionization parameter is expected to be dependent on the 
metallicity because stellar atmospheres of massive O stars become 
cooler with incre asing metallicity as a result of enhanced line and 
wind blanketing (Massev et al. 2005), decreasing consequently the 
ionization parameter. Moreover, when stellar atmosphere abun- 
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Figure 9. Relation between the ionization parameter U and the 
[S n]/[S m],[S iiJ/Ha and [Ora]/[S n] line ratios. The colorful solid lines rep- 
resent results for different metallicities, such as in Fig. [2] and the black solid 
line represents the linear fitting of the average of these results. The dotted 
line represents the relation proposed by Diaz et al. ( 199lT). 



dance is higher, this scatters the photons emitted from the pho- 
tosphere more efficiently, causing a greater conversion efficiency 
from luminous energy flux to mechanical energy flux in the stellar 
wind b ase region, which al so leads to a diminution of U in the H n 
region ( Do pita et al. l2006h . A decrease of U with the i ncrease of Z 
was fo und for H n galaxies bv lNagaoetai]d2006l) and lMaier et al] 
d2006h and for disk H n regions bv lBresolin et al.l ( fl999i) . However, 
our results indicate no systematic dependency of the ionization pa- 
rameter on the metallicity for the sample of objects considered, al- 
though using our grid of models it is impossible to check what fac- 
tor (e.g. stellar effective temperature, geometrical factors, aging) is 
responsible for this behavior. This resu lt is in consonance wi th the 
ones found by iDors & Copettil d2005t) . iGarnett et all d 19971) . and 
Kennicutt & Garnett 119961) . who used the [S n]/[S in] in order to 
estimate the ionization parameter in Hn regions located in spiral 
disks. To analyze whether our result is not an artefact of the meth- 
ods used in this paper, in Fig. [8] we plotted the [S n]/[S m] line ratio 
intensity and the O/H abundances obtained by T e -method of our 
sample. Again, we can note no systematic behavior of U with O/H. 

As noted in Fig. |6] the [Om]/H/J vs. [Nn]/[On] and 
[Oiii]/[Oii] vs. [Nii]/[Oii] indicate very high U values (log U > 
-2) for some objects, which are not predicted by the detailed model 
and by the ([Om]/H/?)/([N n]/Hor) vs. [S n]/[S m] diagram. This oc- 
curs because photoionization models undeipredict the [Om]/Hy6 
(see Fig. [2j, mainly for objects with low metallicity. Because these 
line ratios are age dependent, it is possible that models assuming 
a harder spectral energy distribution would resolve this problem. 
Thus, we ran a grid of photoionization models (not shown) us- 
ing as ionizing source a cluster of 1 Myr and computed U and Z 
values. We found that, for the diagram with [Om]/[On], this new 
grid yields values of U in consonance with the ones from other 
diagrams. However, for the diagram with [Om]/HyS, although a 
better match between the models and the observational data was 
obtained, the U estimations continue to be overestimated in rela- 
tion to the ones via other methods. The O/H estimates practically 



did not change in these cases. A similar problem was also pointed 
out by Stasihska & Izotov (2003), who using a sequence of pho- 
toionization models to reproduce observational data of H n galax- 
ies, which found that the models undeipredict [O ra]/H/3 for objects 
with 7.4 < log(0/H) < 8.0. These authors invoked several mech- 
anisms to explain this discrepancy, such as secondary ionization 
by X-ray s and shocks but a defi nitive conclusion was not reached. 
Recently, iLevesque et al] d2010h compared observational data of a 
large sample of star-forming galaxies with grid of photoionization 
models, such as the ones presented in this paper, but considering 
as ionizing source stellar cluster with different ages and formed 
by instantaneous and continuous star formation. Similar to our re- 
sults, they found that very young models with 0-1 Myr and formed 
instantaneously reproduce better the observed [O m]/Hy6 than the 
older ones, although better agreement is give n by models adopting 
a con tinuous star formation history (see also IPerez-Montero et aU 
2010). These authors did not compare U estimates obtained by dif- 
ferent emission line ratio, but it is probable that the discrepancy 
found by us is maintained even using t heir models with continuous 
star formation. ILevesque et al] (12010) showed that the new gen- 
eration of Geneva evolutionary tracks, which include stellar rota- 
tion, produces a SED more prominent in the higher-energy regime 
(A < 230 A) than the one used here, with rotation effects being 
more important at lower metallicity. Thus, it is probable that if 
these SEDs were used in our photoionization models, the predicted 
intensities of the [O m]/H/J would be larger and U estimates from 
diagnostic diagrams using this line ratio would conciliate with the 
ones from other diagrams. 

As the [S n]/[S m] ratio is weakly dependent on Z, it is useful 
to calibrate it with U, In Fig. [9] we show the U-[S n]/[S m] relation 
predicted by our model s for th e entire range of Z and the relation 
proposed bv lDfaz et al] dl99lh . as well as our results for [S ii]/Hq- 
and [Oiii]/[Oii] versus U. 

A linear fitting of the average of these photoionization model 
results produce 



log U = -1.36 (+0.07) log[SII]/[S III] - 3.09 (±0.05), 
logf/ = -1.66 (±0.06) log[Sn]/Hor- 4.13 (±0.07), 
log U = 1.22 (±0.07) log[OII]/[0 III] - 2.25 (±0.05). 



(1) 

(2) 
(3) 



As can be seen in Fig. [9] the [S iij/Ho- is a good U indicator 
because it shows little variation with the metallicity for Z > 0.2 Z Q 
and uses emission line with near wavelength, being almost inde- 
pendent of the reddening. Our result for U-[S n]/[S m] relation is 
in ver y good agreement with the relation proposed by iDfaz et all 
dl99ll) . 



7 CONCLUSION 

We compared oxygen estimates for a large sample of objects ob- 
tained by direct detection of the electron temperature with those via 
diagnostic diagrams containing strong emission lines predicted by 
photoionization models, as well as from detailed models. Among 
the diagnostic diagrams considered, we found that the ones uti- 
lizing the emission lines [Om]/[On] vs. [Nii]/[Oii], [Om]/H/3 vs. 
[Nii]/[Oii], and [O m]/H/?/[N n]/Ha vs. [S n]/[S in] diagnostic dia- 
grams give O/H values nearest to the T e -method, with differences of 
about 0.04 dex and dispersion of 0.3 dex. Similar results were ob- 
tained using detailed models but with a smaller dispersion, of 0.08 
dex. The origin of the dispersion is probably due to differences be- 
tween the real N/O-O relation of the sample and the one assumed in 
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the models. We did not find any correlation of the ionization param- 
eter with the metallicity for the objects in our sample. We conclude 
that the combination of two line ratio predicted by photoionization 
models, one sensitive to the metallicity and another sensitive to the 
ionization parameter, which does taken into account the physical 
conditions of star-forming regions, gives O/H estimates close to 
the values derived using direct detections of electron temperatures. 
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